Second dip as a signature of ultrahigh energy proton interactions with cosmic 

microwave background radiation 
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We discuss as a new signature for the interaction of extragalactic ultrahigh energy protons with 
cosmic microwave background radiation a spectral feature located at E = 6.3 x 10 19 eV in the form of 
a narrow and shallow dip. It is produced by the interference of e + e~-pair and pion production. We 
show that this dip and in particular its position are almost model-independent. Its observation by 
future ultrahigh energy cosmic ray detectors may give the conclusive confirmation that an observed 
steepening of the spectrum is caused by the Greisen-Zatsepin-Kuzmin effect. 
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Introduction. — The nature and the sources of ultra- 
high energy cosmic rays (UHECRs) are not yet estab- 
lished despite more than 40 years of research. Natural 
candidates as UHECR primaries are extragalactic pro- 
tons from astrophysical sources. In this case, interactions 
of UHE protons with the cosmic microwave background 
(CMB) leave their imprint on the UHECR energy spec- 
trum in the form of the Greisen-Zatsepin-Kuzmin (GZK) 
cutoff Q and a dip [IS, HQ. 

The GZK cutoff is a steepening of the proton spec- 
trum at the energy Eqzk ~ (4 — 5) x I0 19 eV, caused by 
photo-pion production on CMB. This is a very spectac- 
ular effect, but the shape of this steepening is strongly 
model-dependent @, @]. Thus the GZK suppression is 
difficult to distinguish from, e.g., a cut-off due to the 
maximal acceleration energy in a source. The dip is a 
spectral feature produced by p + 7cmb — > p + e + + e~ 
interactions. It is a faint feature, practically not notice- 
able when the spectrum is plotted in the most natural 
way, J P {E) versus E. The dip becomes more pronounced 
in the modification factor [| rj(E) = J P (E) / J™ m (E) , 
where J P {E) is the spectrum calculated with all energy 
losses included, and J p nm (E) is the unmodified spectrum 
calculated with adiabatic energy losses only. The dip is 
clearly seen in the energy-dependence of n(E) and is re- 
liably confirmed by observational data [jj 0] ■ 

In this Letter, we demonstrate the existence of one 
more signature of UHE protons interacting with the 
CMB, which we call the second dip. In many aspects 
it is similar to the first dip. The first dip starts at the 
energy E cq i — 2.3 x I0 18 eV, where pair-production en- 
ergy losses become equal to those due to redshift. The 
second dip starts at the energy £ , cq2 = 6.0 x I0 19 eV, 
where photo-pion energy losses become equal to those 
due to e + e~-pair production. Both features are not seen 
well when the UHECR spectrum is displayed in a nat- 
ural way. While the first dip becomes visible dividing 
the experimental spectrum by the unmodified spectrum 
J p nm (E) oc E~ lg , the second dip appears dividing by the 
smooth universal spectrum (see below). 



Kinetic equation, Fokker- Planck (FP) equation, and 
continuous energy loss ( CEL ) approximation. — We shall 
calculate the diffuse spectrum of UHE extragalactic pro- 
tons assuming a homogeneous source distribution and a 
power-law generation spectrum with spectral index j g . 
In the CEL approximation, the density of UHE protons 
at the present time to can be calculated from the conser- 
vation of the number of protons as 



n p (E, to)dE 



dtQ gcn (E g )dE, 



9 > 



(1) 



where t is the cosmological time and Q gcn oc E g 7s is the 
particle generation rate per unit comoving volume. We 
denote by E g (E,t) the initial energy of a proton gener- 
ated at the cosmological epoch t, if its present (t — to) 
energy is E. The energy evolution E g (E, t) can be easily 
calculated from the known energy losses. The solution 
of Eq. {T]) was explicitly obtained in Refs. 0, @] and for 
a homogeneous distribution of sources it is called uni- 
versal spectrum because it does not depend on the mode 
of propagation, being the same e.g. for rectilinear and 
diffusive propagation Q. The universal spectrum is ob- 
tained in CEL approximation. With higher precision the 
spectrum can be calculated using a kinetic equation, 

dn p 

-P(E,t)n p + / dE'P(E',E,t)n p (E',t) + Q gea (E,t). 

J E 

Here, n p = n p (E,t), H(t) is the Hubble parameter, 
b va ,i-[{E, t) are the energy losses due to pair-production 
treated in the CEL approximation, P(E, t) is the exit 
probability from the energy interval dE due to P7 — > nX 
collisions, and P(E',E,t) is the probability that a pro- 
ton with energy E' produces a proton with energy E in a 
P7 — y nX collision. Introducing x = E/E' and expand- 
ing the regeneration term in Eq. ^ in a Taylor series 
with respect to (I — x), one obtains at order (I — x) 
the CEL equation with the universal spectrum as solu- 



3H(t)n p + A {[H(t)E + b pair (E, t)} n p } (2) 



tion. Including also the (I 



terms, the FP equation 
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emerges as 

3ti 3 

= -SH(t)n p + — {[H(t)E + b tot (E,t)}n p } 



d 2 
dE 2 



[E 2 D(E,t)n p ] +Q gcn {E,t). 



(3) 



Here, btot(E) is the sum of pair-production and photo- 
pion energy losses in the CEL approximation and 
D(E, t) = (SE) 2 /5t is the diffusion coefficient in momen- 
tum space. 




FIG. 1: Upper panel: the trigger function T(E) as func- 
tion of energy. Lower panel: the distortion factor k(E) = 
Jkin{E)/ Juniv(-E) for j g = 2.7 and different values of -Emax. 

The kinetic equation (|2|) allows us a transparent in- 
terpretation of the spectral feature which appears at the 
energy E cq2 = 6.0 x 10 19 eV. At this energy one may 
limit the consideration to the present cosmological epoch 
t k, to. As direct calculations show, the absorption term 
~P(E)n p {E) is then compensated with high accuracy 
by the regeneration term with P(E,E') in Eq. @. At 
E « £ , eq 2, the small CEL term (pair production) breaks 
this compensation, increasing the absorption term, and 
the spectrum acquires a dip. ft is quite narrow because 
photo-pion energy losses increase with energy almost ex- 
ponentially and at E > E cq 2 the pair-production en- 
ergy losses become too small. On the other hand, at 
E < E cq 2 the photo-pion energy losses are too small, the 
spectrum is fully determined by pair-production energy 
losses, while the interference effect disappears. 



Trigger mechanism. — Prior to presenting exact numer- 
ical calculations we shall study semi-quantitatively the 
triggering mechanism, responsible for the second dip. 
One can rearrange the first three terms on the rhs of 
Eq. © into PeffCM = P(E,t)+P cont (E,t) with 



P CO nt(E,t) = 2H{t) 



db(E,t) 



(E,t) 



E 



H(t) 



dE 
d\nn{E,t) 
d\uE 



(4) 



It is the term P con t{E,t) that breaks the above- 
mentioned compensation between absorption and regen- 
eration terms in Eq. ©, triggering thereby the modifi- 
cation of nkin(-E, t). 

It is convenient to introduce the auxiliary trigger func- 
tion T(E) defined at t = to as 



T(E) 



P eS {E)/P cont {E) for E<E C 
P cS (E)/P{E) for E>E C 



(5) 



where E c = 6.1 x 10 19 eV is determined from the con- 
dition P(E) — P cont (E) and is approximately equal to 
E cq2 . The trigger function describes how P c ff(E) is 
changing from P CO nt(E) at E <C E c , where T(E) = 1, 
to P{E) at E > E c , where T{E) = 1 as well. As 
long as T{E) « 1, there is no interference between pair- 
production and pion-production terms, and the ordinary 
solutions are valid. At E — E c , T(E) reaches its max- 
imum and P e g(E c ), being noticeably larger than P(E), 
breaks the compensation between absorption and regen- 
eration terms in Eq. making absorption larger. As a 
result, nkin(-E') decreases around E c . The trigger function 
is plotted in Fig. [T] It reaches its maximum T(E) = 2 
at E = E c « 6.1 x 10 19 eV. As explained above n^i^E) 
must have a local minimum at this energy. The triggering 
mechanism predicts that the position of the dip minimum 
^dip does not depend on j g and £" max and these predic- 
tions are confirmed by our numerical calculations. The 
shape of n^^E) is expected to be similar to the shape 
of the trigger function T(E) and this expectation is also 
confirmed by numerical calculations (see Fig. [1]). 

Numerical solutions. — We next discuss the second dip 
using numerical solutions of the kinetic equation ©. As 
mentioned above, the first dip is distinctly seen in the 
energy dependence of the modification factor rj{E) = 
J p (E)/J™ m (E). Similarly, the second dip is well seen, 
when the spectrum is described by a distortion factor, 
defined as k(E) — Jkin{E)/ J un i v (E), where J U niv{E) is 
the universal spectrum from Eq. |T]). We emphasize 
that the correct prediction for the measured spectrum 
is given by the kinetic equation ©, while the univer- 
sal spectrum, used as a reference spectrum, is obtained 
in the CEL approximation and as such does not include 
the interference between pair and pion production. The 
calculated distortion factor is shown in the lower panel 
of Fig. Q] for 7 g = 2.7 and four values of E max . The 
second dip is clearly seen. Its minimum is given by 
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^2dip = 6.3 x 10 19 eV in good agreement with the predic- 
tion of the triggering mechanism. The width of the dip 
also agrees well with that of the trigger function T(E). 
The independence of the spectral shape of the dip from 
the numerical value of E max , seen in Fig. [T] is another 
prediction of the triggering mechanism. 

The distortion factor k does not return to unity af- 
ter the second dip, but continues to grow for E ^> E c . 
This deviation from unity is explained by fluctuations 
in photo-pion interactions. For E — > £ max , the ratio 



K = Jkin(E)/J C ont(E) 

becomes stationary, 



oo: At these energies Eq. @ 



-P(E)n p {E)- 



dE'P(E', E)n p {E') + Q gcn (E) = . 



When E approaches E max , the regeneration 
term disappears, and one obtains n\ti n (E max ) = 
Q g en(E max )/P(E max ), remaining finite at E max . Using 
the CEL approximation, the equation reads 



8/dE [b tot (E)n p (E)} + Q gen (E) = 0. 



(6) 



For E -> £ max , 

T^cont 

(E) cx (£ max - E)/E max -> and 
hence n(E max ) = n kin (E max ) / n cont (E max ) -> oo. The 
explosive behavior of the distortion factor for E — ► E max 
reflects the different limiting density of n p (E) in the ki- 
netic and CEL equations. This result has a clear physi- 
cal meaning. A particle has a finite probability to travel 
a finite distance without losing energy. While the ki- 
netic equation describes correctly this effect, in the CEL 
approximation particles lose energy for any distance tra- 
versed, however small it may be. This influences the ratio 
n^in(E) / n con t(E) at all energies close enough to E max as 
can be seen in Fig. Q] 

We have obtained one more proof for the energy-losses 
interference as the origin of the second dip. For this we 
have calculated the distortion factor in a toy model in 
which pair production and adiabatic energy losses were 
switched off. Then the interference term must disappear 
together with the second dip. The numerical calculations 
have confirmed this prediction for different 7 S and E max . 

In Fig. [21 the distortion factor is shown for different 
values of j g . One may observe the universality of the sec- 
ond dip with respect to variations of the spectral index, 
as expected from the triggering mechanism. We have 
performed these calculations using the FP equation. 

Figures 1 and 2 show that the second dip is not sensi- 
tive to the exact values of 7 9 and -E max . This implies also 
that a distribution of j g and -E max values does not change 
the shape and position of the dip. Moreover, the cosmo- 
logical evolution of sources is negligible at the energy of 
the second dip. The presence of nuclei primaries affects 
the second dip only for an extreme assumption about the 
fraction of nuclei. Light nuclei are photo-disintegrated 
at this energy and only the heaviest nuclei like Al and 
Fe survive. Their fraction at the production should be 
higher than 20% in order to hide the second dip 12 1. 
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FIG. 2: The distortion factor « as function of E for E max = 
1 x 10 23 eV and different values of the spectral index j g . 



As our final test for the second dip we calculate the 
spectrum with a Monte Carlo (MC) simulation. The re- 
sult of a MC simulation must coincide with the solution 
of the kinetic equation, if all relevant parameters of the 
problem are identical and when the number of MC runs 
tends to infinity. In our case, equal conditions means 
a homogeneous source distribution, the same generation 
spectra and E max as well as identical pj interactions. We 
run the Monte Carlo simulation as described in Ref. Q 
using SOPHIA [l(| for the photo-pion interactions, while 
in the kinetic equation approach the calculations from 
Ref. @ were used. As long as only average energy losses 
are concerned, the results of both works coincide very 
well (see Q for a comparison). However, already small 
differences in the modeling of (differential) cross sections 
of order of a few percent can result in sizable variations 
of the distortion factor k. Numerical errors in the calcu- 
lations are another source of possible discrepancies. In 
Fig. [3] (top), we compare the distortion factors calcu- 
lated with the kinetic equation, FP equation and MC 
methods for an homogeneous source density. The nar- 
row second dip with minimum at i?2dip = 6.3 x 10 19 eV 
is present in all calculations with small differences in 
shapes. The points from MC simulation are connected 
by straight lines, which helps to see the statistical uncer- 
tainties present especially at high energies. 

We have also performed MC calculations for a dis- 
crete distribution of the sources using the values n s = 
1CP 5 Mpc~ 3 inspired by small-angle clustering and the 
very low density n s — 10 -7 Mpc -3 , both shown in the 
bottom panel of Fig. [3] In the first case, the dip agrees 
well with those shown in the upper panel, while in the 
latter case the large distance ~ 200 Mpc to the nearest 
sources results in an early, very steep GZK cutoff that 
covers up the second dip. 

We have compared these calculations with the AGASA 
data [Tll| - The experimental distortion factors are ob- 
tained dividing the observed flux by the universal flux 
and normalizing the distortion factor at low energies to 
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FIG. 3: Top: The distortion factor calculated for -E max = 
1 x 10 23 eV and j g = 2.7 using the kinetic equation (curve 
'kin'), the Fokker-Planck equation (curve 'FP') and Monte 
Carlo calculations (curve 'MC') for a homogeneous source dis- 
tribution together the AGASA data with error bars reduced 
by a factor 3. Bottom: The distortion factor for a discrete 
source distribution, n s — 10 -5 Mpc -3 , inspired by small-scale 
clustering, and n s — 1CP 7 Mpc" 3 . 



K = 1. We have diminished the true AGASA error bars 
by a factor three to give an impression of the potential of 
the Pierre Auger Observatory (PAO) to observe the sec- 
ond dip. This factor corresponds to a factor 10 improve- 
ment in statistics of PAO compared to AGASA. Even if 
the two data points at 4 x 10 19 eV and 6 x 10 19 eV would 
lie exactly on the predicted dip (this is quite possible for 
the true AGASA error bars), the large error bars in the 
PAO data will prevent a reliable conclusion on the pres- 
ence of the second dip. The second dip is expected to 
be seen in the future JEM-EUSO space experiment [l3j |. 
which will have a 100 times higher statistics than Auger 
(see bottom panel of Fig. [3]). However, this expectation 
depends critically on the final energy threshold of this 
experiment, which is currently estimated as 5 x 10 19 eV 
but is planned to be lowered (l3| . The second dip may 
be used as energy calibrator for this experiment, but ac- 
curate MC detector simulations are needed for this con- 
clusion. 

Conclusions. — We have found a new signature of the 
interactions of extragalactic UHE protons with the CMB 
radiation — the second dip. It is explained by the inter- 
ference of pair and photo-pion production and has the 
shape of a narrow and shallow dip. The second dip is 
not seen if the admixture of nuclei heavier than Al in the 
generation flux is larger than 20% and if the space den- 
sity of the sources is extremely small. The observation of 
the second dip is challenging for present experiments such 
as the PAO and requires future high-statistic cosmic ray 
detectors. Combined with the steepening of the UHECR 
spectrum, its observation provides an unambiguous sig- 
nature of the GZK effect. 
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